Lassoing Line Outages in the Smart Power Grid 



Hao Zhu and Georgios B. Giannakis 
University of Minnesota, Dept. of ECE, 200 Union Street SE, Minneapolis, MN 55455, USA 



o : 



00 

o 



> 

00 

rn 
o 



X 



Abstract — Fast and accurate unveiling of power line outages 
is of paramount importance not only for preventing faults that 
may lead to blackouts, but also for routine monitoring and control 
tasks of the smart grid, including state estimation and optimal 
power flow. Existing approaches are either challenged by the 
combinatorial complexity issues involved, and are thus limited 
to identifying single- and double-line outages; or, they invoke 
less pragmatic assumptions such as conditionally independent 
phasor angle measurements available across the grid. Using only 
a subset of voltage phasor angle data, the present paper develops 
a near real-time algorithm for identifying multiple line outages 
at the affordable complexity of solving a quadratic program via 
block coordinate descent iterations. The novel approach relies 
on reformulating the DC linear power flow model as a sparse 
overcomplete expansion, and leveraging contemporary advances 
in compressive sampling and variable selection using the least- 
absolute shrinkage and selection operator (Lasso). Analysis and 
simulated tests on the standard IEEE 118-bus system confirm the 
effectiveness of lassoing line changes in the smart power grid. 

Index Terms — Identification of line outages, compressive sam- 
pling, basis pursuit. Lasso, block coordinate descent. 

I. INTRODUCTION 

It is well appreciated that major blackouts have occurred 
partly due to lack of comprehensive situational awareness 
of the power grid lITSI . This speaks for the importance of 
timely monitoring the status of generators, transformers, and 
transmission lines. Identifying outages and generally changes 
of lines is particularly critical for a number of tasks, including 
state estimation, optimal power flow, real-time contingency 
analysis, and thus security assessment of power systems. 

To appreciate the opportunities and challenges facing the 
line-change identification problem, it is prudent to think of the 
grid as a graph comprising topologically interconnected power 
systems. Phasor measurement units (PMUs) provide voltage 
and power data per local (a.k.a. internal) system in real time. 
Likewise, real-time data are telemetered internally per system 
to offer topology-bearing information on the connectivity 
status of local circuit breakers and switches. On the other 
hand, power flow conservation across interconnected systems 
allows for identifying changes even in external system lines 
- a critical task for comprehensive monitoring. This would 
have been a non-issue if inter-system data were available. 
Unfortunately, the system data exchange (SDX) module of 
the North- American Electric Reliability Corporation (NERC) 
can provide grid-wide interarea basecase topology information 
on an hourly basis |[T6l . but the desiderata is near real-time 
monitoring of transmission lines. In a nutshell, the need arises 
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for each internal system to identify in a computationally 
efficient manner line outages (and generally monitor line 
changes) in its external counterpart relying only on basecase 
topology information and local PMU data. 

Existing approaches typically formulate external line-outage 
identification as a combinatorially complex (integer program- 
ming) problem, which can be computationally tractable only 
for single- or at most double-line outages ||5], ifTTI . ifTZ l. An 
alternative approach adopts a Gauss-Markov graphical model 
of the power network and can cope with multiple outages at 
affordable complexity, but assumes conditionally independent 
phasor angle measurements and requires inter-system PMU 
data to be available in real time across the grid 18]. 

The present paper contributes a computationally efficient 
algorithm for near real-time identification of multiple external 
line outages (and generally changes) using only hourly base- 
case topology information and local voltage phasor angle data 
available by PMUs. It relies on the standard DC linear power 
flow model, but also applies readily to its AC linearized coun- 
terpart per iteration. The novel approach views the topology- 
bearing basecase information as the weighted Laplacian matrix 
of the grid-induced graph (Section |ll}. This leads to an 
overcomplete model of the changes-induced innovation, which 
in turn enables casting line-change identification as a sparse 
vector estimation problem (Section HiH i. Solving the latter 
draws from recent advances in compressive sampling and vari- 
able selection in linear regression problems lH], H, llT3l . and 
leverages the least-absolute shrinkage and selection operator 
(Lasso) to identify line changes at affordable complexity using 
provably convergent block coordinate descent (BCD) iterations 
(Section IIVI ). Simulated tests corroborate the merits of the 
novel algorithm (Section |V]i, and the paper is wrapped up 
with a concluding summary (Section IVII i. 

Notation: Upper (lower) boldface letters will be used for 
matrices (column vectors); (•)-^ denotes transposition; (•)^ the 
matrix pseudo inverse; 1 the all-one vector; the all-zero 
vector; || ■ \\p the vector p— norm forp > 1; | • | the magnitude 
or cardinality of a set; Si\S2 the relative complement of the 
set 52 in the set Si. 

II. Modeling and Problem Statement 

Consider a power transmission network consisting of N 
buses denoted by the set of nodes Af := {!,..., iV}, and 
L transmission lines represented by the set of edges £ := 
{{m, n)} C AfxAf. Collect all the voltage phasor angles {d„}, 
one per bus n e Af, in the vector 6 G M^; and correspond- 
ingly the injected power variables {Pn} in p. The network 



buses comprise the union TV = Ni IJ Me with Mi H Me = 0, 
where Mi denotes the subset of observable buses in the internal 
system, and Me stands for the unobservable buses of the 
external system. Accordingly, 9 comprises two sub-vectors 
Oi and 6e, collecting phasor angles of voltages measured at 
Mi and Me buses, respectively; and similarly for the p/ and 
TpE partitions of p. Supposing that the network-wide injected- 
power vector p only changes gradually across time, the present 
work aims to unveil line changes, including (possibly multiple) 
line outages anywhere in the network, using data 0/ acquired 
in real time from PMUs. Before explicitly formulating the line- 
change identification problem, it is important to understand 
how the network topology dictates the relationship between p 
and 0, as detailed in the following. 

A. Linear DC Power Flow Model 

Power flow models are useful for determining how injected 
power p flows along all transmission lines. To cope with the 
nonlinear AC power-voltage relationships, the DC power flow 
model offers a simple linear analysis tool, assuming constant 
voltage magnitudes per bus and negligible transmission losses 
ifTSl Sec. 4.1]. Although confined only to linear approximate 
analysis of the actual nonlinear system, the DC power flow 
model turns out to facilitate a variety of power system mon- 
itoring tasks under normal operating conditions, including 
security-constrained contingency analysis, and system state 
estimation; see e.g., fTSl Ch. 11-12]. 

In the linear DC model, the power flowing from bus m to 
bus n through their connecting transmission line is given by 

Pran = -^{O-in - On), V (m, u) e £ (1) 

where Xmn = Xnm denotes the reactance along line {m,n). 
Notice that Pmn can be either positive or negative, depending 
on whether actual power flows from bus m to bus n, or in the 
reverse direction. Flow conservation implies that the power 
injected to bus ji balances all the line flows originating from 
it; that is, 

Pn= J2 P"--^ J2 —{dn~drn) (2) 
rn£A/„ mgAin 

where Mn denotes the set of neighboring buses Unked to bus 
n. Writing (|2]i in vector-matrix form yields 

p = Be (3) 

where the N x N matrix B has its (m, n)-th entry given by 

{-^/{xrnn) if (m, n) e £ 

^•^<£Af„ l/(2^ni')' if m = n (4) 
0, otherwise. 

Matrix B, relating the voltage-phasor angle 9 to the injected 
power p as in (|3]l, is uniquely determined by the line reactance 
parameters {xmn} of the network {M,£), and topology- 
bearing information provided by the SDX. At this point, it 
is worth noting that each Xmn reactance is present in only 
four entries of B, namely, (m, to), (to, n), (n, to), and {n, n). 



This observation about line's (to, n) presence in B intuitively 
describes the selective effect topology changes due to line 
outages can have on the angle vector 9 in (|3]l. To formulate this 
topological effect concretely, B can be viewed as the weighted 
Laplacian matrix of the graph {M,£). To this end, consider 
the N X L bus-line incidence matrix M, see e.g., ifTTl pg. 56], 
formed by columns {m^j^^j^ of length N. With subscript i 
corresponding to the line (to, n), the column has all its 
entries equal to except the 7Ti-th and n-th, which take the 
value 1 and —1, respectively. The weighted graph Laplacian 
leads to the following representation of the network topology 
matrix as 

i 1 

B = MD^M^ = V —memj (5) 

xe 

1=1 

where the diagonal matrix D^^ has its £-th diagonal entry 1 /xg 
equal to the inverse reactance 1/xmn, if ^ corresponds to the 
Une {m,n). In addition, matrix B is symmetric, and has rank 
A'^ — 1 if the power network is connected, since its null space 
is only spanned by 1; see e.g., ifTTl pg. 469]. Rank deficiency 
of B gives rise to multiple solutions for 9 in (O. To fix this 
ambiguity, one generation bus is typically chosen as reference 
with its phasor angle set to zero, case in which phasor angles of 
all other buses denote their differences relative to the reference 
bus; see e.g., |i8,, pg. 76]. 

B. Unveiling Network Faults of Line Outages 

Given the linear DC power flow model (|3]l, and the afore- 
mentioned reference bus convention, the pre-event phasor 
angles 9 can be uniquely characterized by the injected power 
vector p and the topology-dependent matrix B. Suppose that 
due to changes in the grid, e.g., cascading failures in an early 
stage, several outages occur on lines, collected in the subset 
£ C £. Line outages on the transmission network yield the 
post-event graph iM,£'), witlfl £' := £\£. 

Without loss of generality (wlog) it is assumed that fast 
system dynamics are well damped, and that the system settles 
down to a quasi-stable state following the line outages. The 
possibility of poor system damping can be accounted for by 
low-pass filtering which smooths out system oscillations, as 
detailed in ifTDl . Under these considerations, the linear DC 
model for the quasi-stable post-event network is given by 

p' = p + V = B'9' (6) 

where B' is the post-event weighted Laplacian of {M,£'), 
and noise vector v accounts for the small perturbations be- 
tween pre- and post-event injected power vectors due to, e.g., 
variations in the bus loads. Typically, bus power injections 
can be modeled as independent Gaussian-distributed random 
variables ifTOl . since each bus load comprises a large number 
of individual users connected through their local distribution 
network. Thus, v can be reasonably assumed zero-mean 
Gaussian distributed, with covariance matrix cr^I. 

'As a mnemonic, post-event quantities are denoted with prime, and the 
differences relative to their pre-event counteiparts are denoted with tilde. 



Voltage phasor angles are only available at the subset of 
buses Mi- Thus, unveiling outages amounts to identifying the 
Unes in £, given pre- and post-event phasor angle Oj and 
O'j. It turns out that solving this problem incurs combinatorial 
complexity. Using (|3]l and (|6]l, the ensuing section will pursue 
an overcomplete representation for line outages. But before 
this, a remark is due on comparing the linear DC power flow 
model here with the probabilistic dependence graph in fSl. 

Remark 1 {Comparison with To capture the spatial 

correlation of phasor angle data across the grid, a Gaussian 
Markov random field (GMRF) approach is pursued in IS] 
to model the statistical dependence of phasor angles among 
buses. Conditioned on {Om]meM„ at neighboring buses, the 
GMRF model in |8| assumes that On is Gaussian distributed 
and independent of all other entries in 0. Unfortunately, 
the GMRF model falls short in capturing the actual power 
system operation for the following reasons. First, it assumes 
(conditionally) Gaussian distributed phasor angles. Second, it 
is not clear why conditional independence holds in practice. 
From the flow model ([T]i, the power conservation law, and with 
{Om}m&^f„ given, the phasor angle On is correlated with those 
of their two-hop neighbors. Lastly, as mentioned in SectionH it 
is currently impossible to acquire 0b in real time. In contrast, 
similar to IS), ifTTl and lfT2l . the present paper adopts the linear 
DC flow model, which has been used extensively for various 
power system monitoring tasks. 

III. Overcomplete Representation of Line Outages 

In this section, the line outage identification task is for- 
mulated as a sparse representation problem. To this end, the 
difference B := B — B' denoting the weighted Laplacian for 
the outage lines in £, is written as [cf. ©J 

B = y — m^mf . (7) 

Substituting the pre-event power flow model (|3]l into the post- 
event one in (|6]l leads to 

B6» + v = B'6»' = B6»' B6/'. (8) 

Consider now the phasor angle change vector 6 := 6' — 6, 
and partition it into the available Oi and the unavailable Be, 
each corresponding to buses in J\fj and A/e; and likewise, 
partition columns of B correspondingly to B/ and B^;. With 
these notations, substituting (|7]) into ([8]), yields 

B0 = B/0/ + Bi30£; = B0' + V = ^s^nif + V (9) 

where si := mJO' /xg, Vi 6 £. To identify £' containing a 
given number of line outages, one can exhaustively test over 
all possible line combinations and select the one offering the 
least-squares (LS) fit of (|9]l. Such an approach incurs combi- 
natorial complexity, and has confined existing methods based 
on extensive enumeration of all combinations, to identifying 
single-line outage lH], ifTTl . or at most double-line outages 
after reducing the search space llT2l . 



To bypass this combinatorial complexity, the fresh idea 
here is to consider an overcomplete representation of all line 
outages. With the available y ;= B/0/ from buses in A//, the 
data model in (|9]l can be reduced to a sparse linear regression 
by introducing the L x 1 vector s, whose i-th entry equals si, 
if £ £ £, and otherwise. Thus, ^ can be written as 

y = y Sfmf-BB0£+v = Ms-B£;^£;+v (10) 

eei 

where the incidence matrix M, now viewed as a regression 
matrix, captures all the pre-event transmission lines. Introduc- 
ing M allows the possibility of multiple Une outages. Different 
from (|9|l, the line outage set £ is no longer present in (fTOl i. 
which eliminates the need for enumeration over all possible 
line combinations. 

In essence, the overcomplete representation in (fTOt allows 
one to cast the problem of recovering the subset £ as one of 
estimating the unknown angle vector 9 e and the sparse vector 
s. The key premise is that the number of line outages is a small 
fraction of the total number of lines; i.e., \£\ <C L. This holds 
even for multiple line outages such as those occurring during 
cascading failures, at least in the early stage when only a small 
number of lines start to fail. The same premise is adopted for 
contingency analysis, where usually single- and double-line 
outages are of primary concern; see e.g., llTS] Ch. 11]. Under 
this sparsity constraint, the underlying outage line locations are 
few, meaning that the vector s has only a few nonzero entries. 
In turn, this suggests recovering line outages by regularizing 
the LS criterion with the ^i-norm of s, which promotes sparse 
solutions. Specifically, the pertinent objective becomes 

{s\0^} :=argmin||y-Ms + B£^B||2 + A||s||i (11) 

S,6 E 

where A is a regularization parameter, and its choice will be 
detailed later on. 

Using the overcomplete representation in (ITOt . the line 
outage recovery problem has been formulated as the convex 
quadratic program (QP) in (fill . Hence, the global minimizer 
of (fTTT i can be efficiently obtained using general-purpose 
convex solvers, such as interior-point algorithms ||3] Ch. 11]. 
Instead of these off-the-shelf solvers, the block coordinate 
descent (BCD) method will be adapted in the next section for 
solving ( fTTT l. because it exploits the specific problem structure 
and can find the line outage identification path for a sequence 
of regularization parameters. Compared to previous works 
constrained to at most double-line outages ||5|, ifTTl . IIT2I . the 
proposed formulation overcomes the inherent combinatorial 
complexity with a virtual "outage" at every possible trans- 
mission line. By leveraging the sparse signal representation 
framework, the proposed approach offers the potential to 
achieve high accuracy in recovering faults with more than two 
line outages, while bypassing the computationally prohibitive 
enumeration of combinations for potential line outages. 

Remark 2 {Partial information on v). The internal system 
formed by buses in Mi can also acquire p/ in real time via 
PMU data and state estimates ||5l. In this case, the internal 



system has also available part of the noise vector v/ at the 
internal buses Mi, and can thus exploit this information when 
estimating the vectors s and 9i in ( fTOl i. Compared to ( fTTT i. 
use of this additional information identifies line outages by 
minimizing the £i-norm regularized LS cost for only part of 
the system in ( fTot corresponding to the buses in Me- With 
V/ available, strict equality constraints can be enforced on the 
remaining system corresponding to buses in Mi- As expected, 
this extra information enhances accuracy and improves identi- 
fiability without sacrificing computational efficiency, since the 
problem is still a QP. Given this similarity, the remainder of 
the paper does not consider this extra information on v/. 
Remark 3 (General models and line faults). The DC power 
flow model has broad applicability in practice, since it also 
coincides with the system obtained after linearizing the de- 
coupled nonlinear AC power flow model, see e.g., JS] Sec. 
5.7]. Thus, the proposed overcomplete representation for line 
outages offers the potential to generalize the nonlinear AC 
power flow model via iterative approximation. Additionally, 
the expression of B in terms of the line reactance parameters 
in (Q suggests that the present analysis is also applicable 
to more general line changes. For example, if the ^-th line 
experiences a sudden change in its reactance, then this can 
be also captured by B as differences in the diagonal weights 
for the Laplacian matrices B and B'. In this sense, the line 
outage is just a special type of line change, since equivalently 
the line reactance goes to infinity, and its corresponding weight 
on matrix B' becomes zero. Therefore, the idea of a sparse 
overcomplete representation s of line outages readily extends 
to the nonlinear power flow model, and can be used to unveil 
general line parameter changes. 

TV. Line-Outage Identification Path via BCD 

Although the QP problem in (fTTT i is efficiently solvable 
using general-purpose convex optimization routines, tuning up 
the regularization parameter A is an issue to be addressed. In 
this section, first the identification path of in ( fTTT i entailing 
line-outage vector estimates with variable sparsity levels is 
obtained using the block coordinate descent (BCD) method 
lfT4l . Subsequently, options are offered for selecting A. The 
so-termed regularization path of Lasso-based solutions as a 
function of A has been studied extensively in the context of 
generalized linear model fitting; see e.g., ||6]. This approach 
is advocated especially for sparse unknown vectors of high 
dimension. Consider for now (fTTt with a fixed regularization 
parameter A. The BCD method optimizes the regularized LS 
cost in (fTTT i by cyclically minimizing over the coordinates, 
namely, the vector Be, and all the scalar entries of s. It 
yields successive estimates of one coordinate, while the rest 
are fixed. With the iterate s(i — 1) for s available at the [i — l)- 
st iteration, V i > 0, the iterate OeH) of 9e can be estimated 
by solving the LS problem as 

2 



: arg mm 



y - Ms{i - 1) - 
- (Bi;)^ [y - Ms(* - 1)] 



'BeOe 



Notice that the pseudo inverse of B^; needs only to be com- 
puted once for all iterations. Once 0b (i) becomes available, 
it remains to update the scalar entries of s(i) for iteration i. 
Suppose that the ^-th entry se{i) is to be found. Precursor 
entries {si{i), . . . , se^i{i)} have been already obtained in the 
i-th iteration along with obtained in closed form as in 

(fT2T l. and postcursor entries {se+i{i — 1), . . . , SL(j — 1)} are 
also available from the previous {i — l)-st iteration. Thus, the 
effect of these given entries can be removed from y by forming 

e-i L 
ei{i) :=y + BEeEii) - ^mjSj{i) ^ mjSj{i-l). 

j=i j=e+i 

(13) 

Using ( fT3T l. the vector optimization problem in (fTTt reduces 
to the following scalar one with se{i) as unknown: se{i) = 
argminsj||ef(i) — mfs^||2 + A|s^|]. This is known to be the 
scalar Lasso problem, which admits a closed-form solution 
expressed via a soft thresholding operator as (see e.g., |16]) 



Si{i) = sign{mj ei{i)) 



m 



A 



|m,||i 



2||m,| 



1,. 



..,L 

(14) 
X, if 



(12) 



where sign(-) denotes the sign operator, and [x\ + 
X > 0, and zero otherwise. 

Cycling through the closed forms (fT2T l- (fT4T l explains why 
BCD here is faster than, and thus preferable over general- 
purpose convex optimization solvers of (fTTt . The LS regres- 
sion in (fT2T l involves only matrix multiplication, while the soft 
thresholding operation in (fT4T l is even faster since the vector 
has only two non-zero entries. Further, the BCD iteration 

is provably convergent to the global optimum {s^, 6^} in (fTTt . 
as asserted in the following proposition. 

Proposition 1: (Convergence of BCD) Given the parameter A 
and arbitrary initialization, the iterates {OE{i), s(*)} given by 
(I12l l and ( 114b converge monotonically to the global optimum 
{s^,^^} of the line outage identification problem in dl 11 1. 

Proof: The argument relies on the basic convergence 
result in [ 14 1. It is easy to check that the two summand terms in 
the cost of (fTTt satisfy the assumptions (B1)-(B3) and (C2) in 
|fT4l . Convergence of the iterates s(i)} to a stationary 

point thus follows by appealing to both Iil4i Thm. 5.1] and 
the regularity property of the cost function lfr4T Lemma 3.1]. 
Convexity of the problem (fTTT i further ensures that every 
stationary point is indeed the global optimum. Monotonicity of 
the convergence follows simply because the cost per iteration 
is strictly non-increasing. ■ 
The identification path is further obtained by applying the 
BCD method for a monotonically decreasing sequence of A 
values. Larger values of A in ( fT4T l force more entries of s{i) to 
be shrunk to zero. Hence, if a large enough parameter A in ( fTTT ) 
is picked, the corresponding will eventually become zero. 
Further, with a decreasing sequence of A values, the optimal 
solution for a large A can be used as a warm start for solving 
( fTTT ) with the second largest A. This way, the line outage 
identification path via BCD exploits both the efficient scalar 



solution in ( fT4l i as well as warm starts, to reduce complexity 
and improve algorithmic stability. The BCD based line-outage 
solution path of ( fTTT i is tabulated as Algorithm [T] 



Algorithm 1 : Input y, M, B^, and a decreasing sequence 
of A values. Output {s^, 0^} in ( fTTT i for each A. 

Initialize with i = —1 and s(— 1) = 0. 

for each A value from the decreasing sequence do 
repeat 

Set i := i + 1 and update the iterate 0b (i) as in (fT2] i. 
for^^ = l,...,Ldo 

Compute the residual ei{i) as in (fT3] l. 

Update the scalar Sf(i) via (fl4l i. 
end for 

until BCD convergence is achieved. 

Save 6^ = 6E{i) and s-^ = s{i) for the current A value. 
Initialize with i = — 1 and the warm start s(— 1) = s(i). 
end for 



A. Tuning the Regularization Parameter 

Algorithm [T] yields the line-outage identification path for 
a sequence of A values. Selection of A is a critical issue 
since a larger A promotes a sparser s^, which translates 
to a smaller number of line outages. With additional prior 
information, existing statistical tests can be adapted to select 
A corresponding to the actual number of line outages. 
Number of line outages is fixed or upper bounded. By 
direct inspection of the regularization paths one can determine 
the value of A, so that the degree of sparsity in equals 
the number of line outages. Thus, the nonzero entries in the 
solution indicate the corresponding outage lines in £. Note 
that the number of Une outages is also assumed known in 
most existing works ||5|, ifTTI . lfT2l . 

If the maximum number of line outages is prescribed, 
the identification path can yield the set £ with all possible 
cardinalities less than the maximum. In order to determine 
the actual cardinality, one can adapt a minimum description 
length (MDL) type test; see e.g., Q- Specifically, the rating 
score for a given number of line outages, becomes the sum 
of the LS reconstruction error achievable with the 8 of the 
specified cardinality, plus a penalty term linearly growing with 
the number of line outages. The "best" number of line outages, 
and thus £, is then selected as the one yielding the minimum 
rating score. 

Variance of injected power noise is known. If the variance cr^ 
of the entries in v is known, one can proceed as follows. For 
each cardinality of the set £, find the subset of buses which 
are not related to those outage lines. The effect of Ms at the 
corresponding entries of y in ( fTOl i then becomes zero. The LS 
reconstruction error from these entries of y can be calculated 
to obtain the sample variance of v, as an estimate of cr^. This 
test selects the set £ with the minimum score, in terms of the 
sample variance deviation relative to the given a^. 
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Fig. 1. Line-change identification patli for witfi a decreasing sequence 
of A values. 



V. Numerical Tests 

The proposed sparsity-promoting fault diagnosis algorithm 
is tested in this section using the IEEE 118-bus benchmark 
system Ig). As in the network is partitioned into internal 
and external systems Mi := {1 - 45, 113, 114, 115, 117} and 
Me := {46 - 112, 116, 118}, respectively. Also, the standard 
deviation for the noise v is set to be 5% of one entry of the 
pre-event injected power vector p. 

The software toolbox MATPOWER QH is used to gen- 
erate the phasor angle measurements as well as the per- 
tinent power flows. The set of outaged branches £ 
{(42, 49), (38, 65), (69, 75)} contains the lines indexed by 
£ = 66,95,115, and is excluded from solving the post- 
event power flow. Notice that all these three outaged lines 
connect to buses in the unobservable external system Me, and 
especially the line (69, 75) links two buses with unobservable 
measurements. 

Algorithm [U is applied to obtain the line outage identifica- 
tion path for with an exponentially decreasing sequence 
of 20 A values, as depicted in Fig. [T] Clearly, a larger A 
yields fewer nonzero entries in s^, and thus leads to less line 
outages. In order to find the best number of line outages, the 
tests based both the MDL criterion and the sample variance 
deviation are implemented for \£'\ = 1 to 5 line outages 
along the identification path. The MDL test scores follow from 
the one for general linear regression problems in [j7j, while 
the variance deviation scores are computed as the absolute 
differences between the sample and actual variances. The 
test scores for different numbers of line outages are scaled 
relative to the maximum ones in either test, and are plotted 
in Fig. 121 The test results demonstrate that the minimum 
scores are achieved at the actual number of line outages 
\£'\ ~ 3. Fig. [3] illustrates the absolution value of entries in s^, 
corresponding to a A value that yields three line outages. All 
entries in Fig. [3] are numerically equal to zero, except for three 
entries, coinciding exactly with the lines in £'. All simulation 
results manifest the effectiveness of the proposed algorithm, in 




Fig. 2. Scaled scores versus a variable number of line outages for both MDL 
and sample vaiiance tests. 
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Fig. 3. Absolute value of versus the line index £, for the A value yielding 



three line outages. 



terms of unveiling multiple line outages and inferring topology 
changes at unobservable external systems. 

VI. Conclusions and Current Research 

A novel line-change identification algorithm was devel- 
oped in this paper for power system monitoring, through an 
overcomplete representation capturing the sparse innovation 
manifested by line changes on phasor angle measurements. 
This new formulation allows identification of multiple (even 
cascaded) transmission line outages, by regularizing an LS 
criterion with the £i-norm of the sparse vector comprising 
the overcomplete representation parameters. As a result, the 
combinatorially complex line outage identification problem is 
converted to a computationally tractable convex QP. Thanks 
to its simple closed-form updates, the resultant BCD-based 
near real-time solver can efficiently compute the entire line- 
change identification path for any prescribed degree of spar- 
sity. Practical statistical tests were also adapted to identify the 
actual number of line outages. Numerical tests demonstrated 
the merits of the proposed scheme in unveiling multiple line 
outages occurring even at the unobservable external system. 



The proposed algorithm is currently evaluated on more 
complicated power testing systems, and various cases of line 
outages. Further enhancements to its accuracy are pursued 
by incorporating information on the injected power at the 
internal system, and also by extending it to cope with network 
faults other than line outages, e.g., changes in transmission 
line parameters. At the same time, the novel overcomplete 
representation for line outages is tested in the context of the 
nonlinear AC power flow model, which has wider applicability 
in practice. 
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